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The magnetic structure and phase diagram of the layered ferromagnetic compound Fe 3 GeTe 2 
has been investigated by a combination of synthesis, x-ray and neutron diffraction, high resolution 
microscopy, and magnetization measurements. Single crystals were synthesized by self-flux reac¬ 
tions, and single crystal neutron diffraction finds ferromagnetic order with moments of l.ll(5)/is/Fe 
aligned along the c-axis at 4K. These flux-grown crystals have a lower Curie temperature Tq «150K 
compared to crystals previously grown by vapor transport (Tc=220K). The difference is a reduced 
Fe content in the flux grown crystals, as illustrated by the behavior observed in a series of polycrys¬ 
talline samples. As Fe-content decreases, so does the Curie temperature, magnetic anisotropy, and 
net magnetization. In addition, Hall effect and thermoelectric measurements on flux-grown crystals 
suggest multiple carrier types contribute to electrical transport in Fe 3 _a;GeTe 2 and structurally- 
similar Ni 3 _ a; GeTe 2 . 


I. INTRODUCTION 

The extensive research on graphene and ultra-thin 
transition metal dichalcogenides (TMD) is naturally 
progressing into studies of van der Waals (VDW) 
bonded heterostructures and application-oriented 
configurations. 1-3 VDW heterostructures are produced 
from layer-by-layer stacking of monolayer (or few-layer) 
materials, which can be obtained from exfoliation of bulk 
sources. 4 These engineered structures will inevitably 
lead to unique properties and potentially new technolo¬ 
gies, in part due to the ability to combine materials 
with complementary functionalities. This was recently 
demonstrated with the realization of highly-efficient, pi¬ 
cosecond photoresponse in a heterostructure of graphene 
(fast response) and WSe 2 (high efficiency). 5 The con¬ 
struction of VDW heterostructures using topological 
insulators (TI) will certainly provide ample phase space 
for exploring novel states of matter and theoretical 
predictions. Currently, there is interest in controlling 
the dispersion and carrier type of the Dirac surface 
states of a TI using van der Waals heterostructures that 
allow separate tuning of bulk and surface states. 6 

VDW-bonded ferromagnets are of interest as building 
blocks for heterostructures designed for use in spin-based 
information technologies, either for the direct exploita¬ 
tion of their magnetic properties or via magnetic prox¬ 
imity effects. The latter permits the use of nominally 
non-magnetic materials in spintronics, and is being pur¬ 
sued using EuO/graphene heterostructures. 7,8 Similarly, 
skyrmions are of interest from a fundamental perspec¬ 
tive and for their ability to potentially enable low-power 
spintronics, and these spin states are stabilized by a re¬ 
duction from three to two dimensions and by the presence 
of Rashba spin-orbit coupling. 9,10 In general, spin-orbit 
coupling within heterostructures should yield interesting 
spin structures and magnetoelectric transport. 

Of the ‘next-generation’ VDW materials, Crl 3 and 


CrSiTe 3 -type compounds have recently been identified as 
promising systems with the potential for long-range mag¬ 
netism in monolayers. 11-15 These materials are ferromag¬ 
netic at approximately 61K and 33 K, respectively, 14,16 
which is somewhat low for spintronic applications. Inter¬ 
estingly, magnetic odering temperatures have been pre¬ 
dicted to increase when monolayers are constructed, 12 
and initial experimental results seem to confirm this be¬ 
havior in CrSiTe 3 . 1( Considering this, the VDW-bonded 
compound FesGeTe 2 may be of particular interest be¬ 
cause the bulk is ferromagnetic near 230 K. 18 

Fe 3 GeTe 2 contains FesGe slabs separated by VDW- 
bonded Te layers, as shown in Fig. 1(a). The FesGe 
slabs contain Fe(l)-Fe(l) pairs across a hexagonal net¬ 
work built by Fe(2)-Ge, and this is structurally similar 
to the more three-dimensional Fe 2 _a,Ge compounds. 19 
Fe 3 GeTe 2 is an itinerant ferromagnet, with Curie tem¬ 
peratures of 220 and 230 K reported, and an estimated 
spontaneous magnetization of 1.6/r_e/Fe at OK. 1,18 Pre¬ 
vious anisotropic magnetization measurements on crys¬ 
tals grown using chemical vapor transport suggested that 
the c-axis is the easy axis, and an anisotropy field of 
at least 5T was demonstrated at 10 K. 1 The Fe(2) po¬ 
sition was reported to have a small concentration of va¬ 
cancies (17%), but chemical characterization suggested 
the composition is Fe 3 GeTe 2 . 18 Interestingly, Fe vacan¬ 
cies in Fe 2 _ s Ge are also concentrated on the chemically 
similar Fe(2) position and a wide, complex phase-width 
is observed. 19 A modest phase-width is thus expected in 
FesGeTe 2 based on the published crystallographic data, 
its intermetallic nature, and structural similarities with 
Fe 2 _ 2 ;Ge compounds. Upon completion of this work, the 
existence of a phase-width was independently confirmed, 
though the influence on the structure and magnetic prop¬ 
erties was not reported. 21 

We report the growth of single crystals via a molten 
flux technique, together with a study of polycrystalline 
samples that confirms a phase width and corresponding 
response in the magnetic properties. We find that the lat- 
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tice parameters, Curie temperature and saturation mag¬ 
netization vary smoothly with Fe concentration, though 
the a and c lattice parameters trend oppositely. Reduced 
Fe content in the flux grown crystals is found to be re¬ 
sponsible for their lower Tq when compared to prior re¬ 
ports. Additionally, we find that Fe 3 _ x GeTe 2 itinerant 
ferromagnets possess multi-carrier electronic transport. 


II. EXPERIMENTAL DETAILS 

Single crystals were grown from Fe-Ge-Te ‘self’-fluxes, 
and a melt composition of Fe 2 GeTe 4 produced the largest 
crystals of those investigated. Crystals were also suc¬ 
cessfully grown from the compositions FeGe 2 Te 4 and 
FeGeTe 2 using similar heating procedures, though these 
crystals were typically smaller and had slightly lower 
Curie temperatures. High purity elements from Alfa Ae- 
sar (Fe 99.98%, Ge 99.9999%,'and Te 99.9999%) were 
combined in AI 2 O 3 crucibles and sealed in evacuated 
quartz ampoules. A crucible filled with quartz wool was 
placed on top of the growth crucible to catch the excess 
flux during centrifugation. The melt was homogenized 
at 950°C for approximately 12 h, then cooled slowly to 
675°C, at which temperature the ampoules were removed 
from the furnace and placed in a centrifuge to expel the 
excess flux. A variety of cooling rates were found to pro¬ 
duce crystals, and in this paper we report data from crys¬ 
tals obtained after cooling at 1 and 3°/h (the properties 
were observed to be equivalent). We also synthesized 
crystals of Ni 3 GeTe 2 to provide a non-magnetic reference 



FIG. 1. (color online), (a) Image of crystal structure with 
atomic positions labeled; Fe(l)-Fe(l) bonds pierce the center 
of each hexagon formed by Fe(2)-Ge (space group PQ^/mmc 
- no. 194). (b) X-ray diffraction data for an as-grown facet 
of Fe 3 GeTe 2 showing the c-axis normal with exclusively 001 
diffraction peaks observed; a Le Bail fit (solid line) yielded the 
c-axis lattice parameter shown in the image, (c) A picture of 
a large single crystal. 


material during our investigation of transport properties. 
Crystals of Ni 3 GeTe 2 were grown from a melt of NiGeTe 2 
cooled at 3°/h, with the flux removed at 660° C. The 
crystals of Ni 3 GeTe 2 were generally smaller than those 
of Fe 3 GeTe 2 - Upon characterization, both were found to 
be transition metal deficient. 

Polycrystalline samples with nominal Fe concentra¬ 
tions between Fe 3 .ioGeTe 2 and Fe 2 . 6 oGeTe 2 were pre¬ 
pared by grinding Fe powder, Ge powder, and Te shot in 
a He filled glove box. The mixture was transferred, un¬ 
der He, to a vacuum line where the quartz ampoule was 
sealed under vacuum. The samples were heated at 675°C 
for approximately 10 d. The as-reacted Fe 3 _ 2 ;GeTe 2 sam¬ 
ples were found to be either slightly sintered and dull 
black for small x , or heavily sintered with visible grain 
growth/crystallization for large x. At this point, the sam¬ 
ples were analyzed with x-ray diffraction and magnetiza¬ 
tion measurements. Small pellets were then fired briefly 
at 600° C to facilitate isothermal magnetization and en¬ 
ergy dispersive spectroscopy (EDS) measurements. EDS 
measurements were performed using a Bruker Quantax 
70 EDS system on a Hitachi TM-3000 microscope. 

Crystals from the flux-grown reactions are much larger 
than those needed for single crystal x-ray diffraction. 
Due to the ease with which the crystals cleave and de¬ 
form, the process of cutting or crushing the crystals in¬ 
duces significant damage. Therefore, a small crystal from 
the polycrystalline reaction with nominal composition 
Fe 2 . 75 GeTe 2 was selected for single crystal x-ray diffrac¬ 
tion. Data were collected at 173 K on a Bruker SMART 
APEX CCD, using Mo -Ka radiation (A = 0.71073 A, 
graphite monochromator). For refinement of the crystal 
structure, absorption corrections were applied with SAD- 
ABS and SHELXL-97 was used to refine the data, and 
the atomic coordinates were standardized with Structure 
Tidy within PLATON. 22-24 The refinement utilized 138 
unique reflections from 1885 reflections and 12 refine¬ 
ment parameters. The Goodness of Fit was 1.303 while 
Rint = 0.0280. Powder x-ray diffraction data were col¬ 
lected using a PANalytical X’Pert Pro MPD with a Cu 
K a 1 (A=1.5406A) incident beam monochromator, and 
Le Bail and Rietveld refinements were performed using 
FullProf. 25 

Magnetization measurements were performed in a 
Quantum Design Magnetic Property Measurement Sys¬ 
tem, as well as with the AC Magnetic Susceptibility Op¬ 
tion on a Quantum Design Physical Property Measure¬ 
ment System (PPMS). A PPMS was also used to mea¬ 
sure the Seebeck coefficient, thermal conductivity and 
electrical resistivity using the Thermal Transport Op¬ 
tion. For this measurement, gold-coated copper leads 
were attached to the crystals using H20E Epo-Tek sil¬ 
ver epoxy. Hall effect measurements and magnetoresis¬ 
tance measurements were performed using a standard 
four-point configuration with Pt wires attached via sil¬ 
ver paint (DuPont 4929N). The Hall resistance ph was 
obtained via the odd-in -H part of the transverse resis¬ 
tance pn = ( p xy [H\ — p xy [-H \)/2 with maximum mag- 

























3 


netic fields of magnitude 17=80 kOe applied along the 
c-axis. Magnetoresistance was obtained from the even- 
in -H portion of the longitudinal resistance with fields 
applied along the c-axis. At 2 K and 80 kOe, the magne¬ 
toresistance was less than 2% and we have excluded the 
data from the manuscript. 

To probe the microscopic crystal and magnetic struc¬ 
ture, we performed single crystal neutron diffraction on 
the Four-Circle Diffractometer (HB-3A) at the High Flux 
Isotope Reactor (HFIR), ORNL. A single crystal of ap¬ 
proximate dimensions 4 mm x 4 mm was mounted on an 
A1 rod inside a CCR. Using an incident wavelength of 
1.003A (Si(331) monochromator), measurements were 
performed between 4K and room temperature. A large 
number of reflections were collected at 220 K and 4 K to 
determine the nuclear and magnetic structures, respec¬ 
tively. Neutron powder diffraction was performed on the 
HB-2A Neutron Powder Diffractometer at HFIR, ORNL. 
A wavelength of 2.41 A was used in all measurements, 
and this was selected with a Ge(113) monochromator; 
the samples were placed in A1 cans. All neutron diffrac¬ 
tion data were refined using the program FullProf. 25 

High resolution scanning transmission electron mi¬ 
croscopy (STEM) was performed using a Nion Ultra- 
STEM200 microscope operating at 200 keV, equipped 
with a Gatan Enfinium spectrometer for the in situ col¬ 
lection of electron energy loss spectra (EELS). Samples 
were examined in both plan view and cross-sectional ori¬ 
entations. Samples were prepared by a combination of 
polishing and Ar + ion milling using a voltage of 2keV. 
Contact with moisture was avoided during sample prepa¬ 
ration. 


III. RESULTS AND DISCUSSION 

The single crystals obtained from self-flux growths 
are thin plates with dimensions reaching greater than 
lcm (see Fig. 1). X-ray diffraction data collected from 
the surface of as-grown facets confirm the expected ori¬ 
entation with [001] normal to the facet, as shown in 
Fig. 1. A Le Bail fit to the diffraction data in Fig. 1 
yielded c=16.3941(2)A. Rietveld refinement of data col¬ 
lected on ground crystals resulted in a=3.9536(7)A and 
c= 16.396(2)A. These values differ sharply from those in 
the literature, where a=3.9910(10)A and c=16.336(3)A 
were reported from room temperature single crystal x-ray 
diffraction. 18 

Ni 3 GeTe 2 crystals were grown from a self-flux to pro¬ 
vide a non-magnetic reference material during the char¬ 
acterization of physical properties. The lattice parame¬ 
ters obtained for our Ni 3 GeTe 2 crystals also differ from 
those in the literature, with refinement of powder diffrac¬ 
tion data yielding a=3.8373(2)A and c=16.048(2)A. 
This compares to a=3.9110(10)A and c=16.022(3)A 
previously reported for crystals with a composition of 
Ni 2 . 95 GeTe 2 obtained from refinement of single crystal 
x-ray diffraction data. 18 Thus, for both Fe 3 GeTe 2 and 
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FIG. 2. (color online). Powder x-ray diffraction data for two 
polycrystalline samples with refined lattice parameters and 
compositions provided. 

Ni 3 GeTe 2 grown from a self-flux, the a lattice param¬ 
eter is smaller and the c lattice parameter is slightly 
larger than the literature reports for crystals obtained 
from nominal 3-1-2 compositions (formed via solid state 
reactions). 

In addition to having different lattice parameters than 
those in the literature, the flux-grown Fe 3 GeTe 2 crys¬ 
tals have lower Curie temperatures than those previ¬ 
ously reported for Fe 3 GeTe 2 - The crystals were grown 
from self-fluxes, which minimize the chance for extrin¬ 
sic doping. We thus speculated that Fe 3 GeTe 2 con¬ 
tains a non-trivial phase width that influences struc¬ 
ture and physical properties. Our EDS results suggest 
Fe and Ge deficiencies may exist, with a composition 
of Fe 2 . 9 i( 3 )Ge 0 g 5 ( 4 )Te 2 .oo( 4 ) obtained. For the Ni-based 
analogue, EDS yielded Ni 2 . 4 o( 4 )Gei.oi( 3 )Te 2 .oo( 3 ) for crys¬ 
tals grown from NiGeTe 2 - 

To investigate the phase width of Fe 3 _a;GeTe 2 , poly¬ 
crystalline samples were produced via solid-state re¬ 
actions with nominal compositions Fe 3 _ :r GeTe 2 with 
—0.1 < x < 0.4. These materials were characterized 
using powder x-ray diffraction, EDS, and magnetization 
measurements. We begin by presenting the results of 
our crystallographic studies, and then return to the mag- 
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netic properties and detailed characterization of our sin¬ 
gle crystals. When necessary for unit conversion, we have 
utilized compositions obtained by Rietveld refinement. 
When discussing single crystals, we present data for crys¬ 
tals grown from Fe 2 GeTe 4 melts, and refer to these by 
the composition obtained from refinement of single crys¬ 
tal neutron diffraction data, Fe 2 ,76Geo.94Te2. 


A. Structure and Composition 

All polycrystalline samples were found to be phase 
pure within the limits of our powder x-ray diffractome¬ 
ter, with the exception being a sample of nominal com¬ 
position Fe 2 . 6 oGeTe 2 - This sample contained GeTe and 
FeTe 2 impurities, and data for this sample have been ex¬ 
cluded from the manuscript. Our results thus suggest 
the Fe-deficient phase boundary is likely reached near 
Fe 2 . 7 GeTe 2 - We have not investigated the phase width 
with regard to Te or Ge. 

The sample of nominal composition Fe 3 .ioGeTe 2 did 
not possess any obvious impurities by x-ray diffraction. 
The refined composition for this sample is Fe 2 . 97 ( 2 )GeTe 2 , 
while large-area (300 /jm diameter) EDS scans yielded 
Fe 3 .o 7 ( 2 )Geo. 92 ( 2 )Te 2 .oo( 2 )- The EDS results are likely in¬ 
fluenced by minor secondary phases that are not easily 
detected by x-ray diffraction. An Fe-rich impurity was 
observed in back scattered electron images, and we were 
not able to isolate large grains of the main phase for EDS 
due to resolution limits. Therefore, we focus on compo¬ 
sitions obtained from Rietveld refinements of the powder 
diffraction data. In the structurally-similar Ni 3 GeTe 2 
compound, occupation of an interlayer Ni(3) position 
at (0,0,0) was detected by single crystal x-ray diffrac¬ 
tion, though occupation of this site was not observed 
in Fe 3 GeTe 2 - 18 We performed refinements including the 
hypothetical Fe(3), and the results suggest that an occu¬ 
pation of 3(1)% may exist for the nominal Fe 3 .ioGeTe 2 . 
This seems reasonable, but the result is dependent on the 
range of data analyzed (in part due to a connection to the 
refinement of sample texture). Based on our single crys¬ 
tal diffraction and electron microscopy, discussed below, 
we feel confident that there is not any significant occu¬ 
pation of this interlayer site in the nominally Fe-deficient 
samples (at least for a;~0.25). We therefore utilized the 
published crystal structure, with two Fe positions, to re¬ 
fine our diffraction data and report compositions based 
on this refinement. 

Figure 2 presents powder x-ray diffraction data and 
Rietveld refinements for samples of nominal composition 
Fe 3 .ooGeTe 2 and Fe 2 . 7 sGeTe 2 , which refined to the com¬ 
positions Fe 2 .g 6 (i)GeTe 2 and Fe 2 . 69 ( 4 )GeTe 2 , respectively. 
The lattice parameters obtained from refinement of the 
data in Fig. 2 confirm that the Fe-deficient sample has 
smaller a and larger c than the Fe-rich sample. This 
trend is shown in Fig. 3 for all polycrystalline samples, 
where the normalized lattice parameters are plotted as a 
function of Fe content. Data obtained from single-crystal 



FIG. 3. Normalized lattice parameters as a function of refined 
Fe content for Fe 3 _ x GeTe 2 samples (powder x-ray diffrac¬ 
tion ‘p" xr d’ at room temperature), including results obtained 
from single crystal x-ray (‘sc-xrd’) diffraction data collected 
at 173 K; a ma x =4.0244(1) Aand c maa; = 16.405(1) A. Vertical 
error bars are smaller than the size of the data markers. 

x-ray diffraction (sc-xrd) are also included to confirm this 
behavior; note that the sc-xrd data were collected at 
173 K where slightly smaller lattice parameters are ex¬ 
pected. The increase in a with increasing Fe may be ex¬ 
pected, with the lattice expanding as Fe(2) vacancies are 
filled. Indeed, a similar trend in a occurs for Fe 2 -a;Ge in 
the range 0.15 < x < 0.7, though c was found to increase 
with increasing Fe content in this more 3D compound. 19 

Given the layered structure, we expect the greatest 
influence of the increased occupation of Fe(2) to be ob¬ 
served within the (Fe 3 Ge) planes. Indeed, the relative 
expansion of a is greater than the relative contraction of 
c (see Fig. 3). We therefore expect the changes in c to 
be driven by the changes in a. That is, the decrease in c 
is due to an expansion of the basal plane that allows Te 
atoms to get closer to the Fe 3 Ge layer, while increased oc¬ 
cupation of Fe(2) also leads to greater Fe(2)-Te bonding. 
Consistent with this, Fe(2)-Te bond distances are small¬ 
est in the Fe-rich samples, as are Fe(l)-Fe(l) bond dis¬ 
tances. As a result of Te being pulled towards the Fe 3 Ge 
layers, the structure collapses along c to maintain the 
van der Waals bonds (the length of which increases with 
increasing Fe). If the hypothetical Fe(3) becomes occu¬ 
pied at high Fe concentrations, its presence could lead 
to increased bonding that would also reduce c. However, 
refinements with Fe(3) reveal a larger Fe(3)-Te bond dis¬ 
tance for higher Fe contents, suggesting that Fe(3) does 
not play a role in the contraction of c. 

We note that the occupancy of Fe(2) is correlated to 
the refinement of texture (preferred orientation) in the 
powder diffraction data. In these samples, the texture 
physically increases with decreasing Fe content. This 
trend can be observed visually, with the Fe-deficient sam¬ 
ples demonstrating more grain growth during the reac¬ 
tion. Not refining the texture results in much larger resid¬ 
uals, and thus the preferred orientation was refined for all 
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samples and care was taken to minimize the texture dur¬ 
ing sample preparation. The main result is not strongly 
influenced by this, though, and the nominal composition 
clearly influences the amount of Fe in the final specimen 
and the trends are consistent. Also, to simplify the re¬ 
finement, we did not allow Ge content to vary and fixed 
an overall displacement parameter. 

TABLE I. Selected data from refinements of single crystal x- 
ray diffraction collected on a crystal from a reaction of nom¬ 
inal composition Fe 2 . 7 sGeTe 2 . Atomic coordinates are Fe(l): 
0,0,2; Fe(2): |,|,A Te: §,§, 2 ; Ge: §,§,§. 

a (A) 3.9421(9) 

c (A) 16.378(5) 

Ri (all data) 0.0461 
wR .2 (all data) 0.1019 
Goodness of fit 1.306 
Fe(l) 2 coord. 0.1721(2) 

Te 2 coord. 0.0900(1) 

The potential existence of interlayer Fe is important 
beyond understanding the modifications to the lattice 
with changing Fe content. As shown below, the total 
Fe content clearly influences the magnetic properties and 
interlayer Fe could play a particularly important role 
in determining the saturation magnetization and coer¬ 
cive field. In regards to the chemistry and structure 
of these materials, the existence of interlayer Fe would 
imply a more three-dimensional material. As such, in¬ 
terlayer Fe would likely impede the production of thin- 
layers by cleaving bulk crystals. Similarly, if an amount 
of interlayer Fe can be controlled, it would likely influ¬ 
ence the anisotropy of the electronic and magnetic prop¬ 
erties. Interestingly, as discussed below, the samples with 
the highest Fe content have the largest anisotropy in the 
magnetic properties. While our powder diffraction re¬ 
sults suggest Fe(3) is not occupied, we felt that additional 
investigation was warranted. 

We investigated the issue of interlayer Fe with single 
crystal x-ray and neutron diffraction, and electron mi¬ 
croscopy. The refinement results from the x-ray diffrac¬ 
tion data are shown in Table I, and we clearly observe 
the partial occupation of Fe(2) with a refined composition 
Fe 2 . 7 i( 2 )GeTe 2 (recall this is for a crystal obtained from 
the Fe 2 . 75 GeTe 2 polycrystalline reaction; Ge vacancies 
were not observed). Electron density was not detected 
at an interlayer Fe(3) position, consistent with earlier 
work for crystals selected from a polycrystalline sample 
of nominal composition Fe 3 GeTe 2 - 18 Our single crystal 
neutron diffraction on flux-grown crystals, discussed in 
more detail below, does not find any strong evidence for 
occupation of the Fe(3) position. 

A large number of vacancies are refined on the Fe(2) 
position from both the neutron and x-ray single crystal 
diffraction data. When filled, these Fe(2) positions are 
bonded to three in-plane Ge atoms. As vacancies are in¬ 
troduced, Fe(2)-Ge bonding decreases (on average) and 
this results in a large in-plane displacement parameter 


TABLE II. Anisotropic displacement parameters for 
Fe 3 _a,GeTe 2 from refinements of single crystal x-ray diffrac¬ 
tion data at T=173K; note Z7n = U 22 = 2 Ui 2 - 


species 

Un 

u 33 

U eq 

occupancy 


A 2 

A 2 

A 2 

fractional 

Fe(l) 

0.0069(12) 

0.014(2) 

0.0091(10) 

1 

Fe(2) 

0.0069(12) 

0.014(2) 

0.0091(10) 

0.71(2) 

Ge 

0.044(2) 

0.015(2) 

0.034(1) 

1 

Te 

0.0081(6) 

0.0152(8) 

0.0105(5) 

1 


for Ge, as shown in Table II. Consistent with this, we 
refine a larger concentration of vacancies on Fe(2) com¬ 
pared to Ref. 18, and our refined Un for Ge is also larger 
than previously reported. However, the previous study 
did report a strong anisotropy for the displacement pa¬ 
rameters of Ge, and a similar U 33 to that shown in Table 
II was reported. 18 In the binary compound Fei. 6 oGe 2 , it 
has been suggested that the Ge atoms actually move off 
of their site and the symmetry is potentially broken. 26 

The data in Table II correspond to refinement of sin¬ 
gle crystal x-ray diffraction data, where only Fe(2) was 
found to be partially occupied. We obtained a fractional 
occupation of 0.71(2) for Fe(2), and this was utilizing a 
constraint to maintain equal displacement parameters for 
Fe(l) and Fe(2). When refined separately, the displace¬ 
ment parameters of Fe(2) are observed to be rather small 
within the basal plane (t/n=0.004(2)). In this case, the 
refined occupancy of Fe(2)=0.69(2) is very similar to that 
obtained from the restrained fit. 

We also performed high resolution scanning transmis¬ 
sion microscopy (STEM) to investigate the local struc¬ 
ture. For imaging, we used the high angle annular dark 
field detector (HAADF), which yields a contrast nearly 
proportional to Z 2 (Z = atomic number) and is therefore 
better suited for imaging vacancies and interstitials, and 
cations with different atomic numbers. 

A representative cross-sectional image with the elec¬ 
tron beam parallel to [130] is shown in Fig. 4(a), and plan 
view images with the electron beam parallel to [001] were 
also collected but are not shown. In general, the crystals 
looked to be of high quality, lacking obvious two or three 
dimensional defects. In the .2-contrast images, the van 
der Waals gaps between layers of Te atoms (brightest 
spots) can be easily identified and intensity is not ob¬ 
served between them, suggesting that interlayer Fe is not 
present in our crystals. In both cross-sectional and plan 
view images, we observe a variation in the intensity of the 
atomic columns containing Fe(2) and Ge, which is consis¬ 
tent with the presence of Fe vacancies. The intensity pro¬ 
files in Fig. 4(b) reveal these variations, with the Fe(2)-Ge 
columns displaying a clear disruption of the pattern ex¬ 
pected for the case of fully occupied Fe(2) sites. Due to 
the difference in Z between Fe (.2=26) and Ge (Z= 32), 
a fully occupied structure would display higher-intensity 
peaks adjacent to lower-intensity peaks in a high-low- 
liigh sequence when viewed down [130]. The red arrows 
in the image highlight Fe and Ge columns with markedly 
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FIG. 4. (color online) HAADF STEM image showing 17- 
contrast. The heaviest columns (Te) show up as the brightest 
spots and Fe columns show up with the dimmest contrast. 
Van der Waals gaps are clearly visible as black stripes. The 
inset shows an overlay of the lattice projection along [130]. 
(b) Profiles of the intensity along the dashed lines shown in 
(a) indicating significant disorder in the Fe(2)-Ge plane, (c) 
EELS integrated intensity maps showing the chemical signa¬ 
ture of the alternating Te, Fe, and Fe-Ge planes. 


lower or higher intensity compared to the average. Based 
on the images, we cannot exclude the possibility of Ge-Fe 
antisite defects or Ge vacancies in the Fe(2)-Ge planes. 
Due to the low vacancy concentration and the more de¬ 
localized nature of the EELS signal as compared to the 
HAADF signal, vacancies are not resolved in the EELS 
compositional maps shown in Fig. 4(c), which show a uni¬ 
form distribution for the integrated intensities of the Te- 
M 4 i 5 , Fe-L 2 ,3 and Ge-L 2j3 edges. 

The images were observed to change upon continued 
exposure to the electron beam. A movement of atoms 
and vacancies could be observed, with atoms eventually 
occupying the interlayer region and vacancies appearing 
to occupy sites besides Fe(2). In order to avoid beam- 
induced hopping and preserve observation of the intrinsic 
structure, images were collected rapidly on unirradiated 
regions by summing 20 frames, each one acquired within 
1 s. EEL spectrum images were acquired using a spac¬ 
ing of 0.87 A and a dwell time of 0.2 s. Beam-induced 
structural effects were also noted for Ni 3 GeTe 2 in Ref. 
18. 

EELS was utilized to look for differences in the oxida¬ 
tion state of the two Fe sites. The ratio of the Fe L 3 j2 
peaks gives an indication regarding the oxidation state, 
though this can vary with details of the bonding environ¬ 
ment. We find an L 3 i2 ratio that varies slightly between 
the Fe(l) and Fe(2) positions, which may suggest the Fe 
atoms are in slightly different oxidations states or carry 
a different moment. 

The lack of evidence for occupation of Fe(3) in the 
STEM and diffraction studies may be due to the low- 
occupation of this position, which is likely to be further 
reduced with the overall decrease of Fe content in our 
crystals due to the flux growth. Examination of Fe-rich 
crystals will offer the best chance for identifying inter- 



FIG. 5. (color online). Rocking curves of the 010 reflection 
from single crystal neutron diffraction on Fe2.76Geo.94Te2 at 
10 K and 157 K. The solid lines are Gaussian fits to the peaks. 
The image on the right shows the magnetic structure obtained 
by refinement of neutron diffraction data collected at 41< (Fe 
= red, Te = green, Ge = blue). 


layer Fe, if it exists, though even then the concentration 
will likely be very small and observation will be challeng¬ 
ing. Even at such low concentrations, though, the pres¬ 
ence of interlayer Fe could strongly influence the physical 
properties. 


B. Neutron Diffraction 

Single crystal neutron diffraction was performed on 
a crystal obtained from the Fe 2 GeTe 4 flux. This flux 
composition was found to produce the largest crystals, 
and thus it was investigated in the most detail (includ¬ 
ing magnetization and transport measurements below). 
The crystal structure was refined using data collected 
at 220 K, and the magnetic structure was obtained from 
data collected at 4K. As discussed below, powder neu¬ 
tron diffraction measurements were also performed to 
verify the magnetic structure and investigate variations 
in the Fe moment(s) as a function of composition. 

Figure 5 shows representative rocking curves obtained 
during the single crystal neutron diffraction data col¬ 
lection. As shown, there is strong temperature depen¬ 
dence that indicates the presence of magnetic order be¬ 
low 157 K. The nuclear refinement is consistent with the 
single crystal and powder x-ray diffraction data. We re¬ 
fine vacancies at the Fe(2) and Ge sites, and do not ob¬ 
serve strong evidence for Fe substituting for Ge. Con¬ 
sistent with the above discussion, we do not detect any 
Fe between the Te-Te layers. The refined composition is 
Fe 2 . 76 ( 4 )Ge 0 . 94 ( 4 )Te 2 with Rf=A.11 and x 2 =0.423. The 
site occupancies were fixed to these values for the lower 
temperature refinements. 

An ordered moment along c of l.ll(5)/xs/Fe is ob¬ 
tained from refinement of the single crystal data at 4K, 
and a schematic of the magnetic structure is shown in 
Fig. 5. Our refinement does not indicate any signifi¬ 
cant difference in the moments on the two Fe positions 
for the Fe-deficient crystal. If we allow both moments 
to refine separately, we obtain 1.07(ll)/Zfl/Fe(l) and 
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FIG. 6. (color online). Integrated intensities of Bragg peaks 
from single crystal neutron diffraction on Fe 2 . 76 Geo. 94 Te 2 - 
The increase in intensity for the 010 peak below Tc originates 
from the magnetic order along c, whereas the intensity of the 
002 peak does not change significantly with temperature. 


0.9(5)/is/Fe(2). We note that the errors on the moment 
values increase appreciably when both are allowed to re¬ 
fine separately, and this gives values consistent with that 
for fixed Fe(l) and Fe(2). Thus, within the limits of 
this data, we have no reason to suspect that the differ¬ 
ent sites carry significantly different moments at this Fe 
concentration. Very recently, based on powder neutron 
diffraction data from Fe 2 .gGeTe 2 , the ratio of moments 
between Fe(l) and Fe(2) was found to be 1.25 at 1.5 K. 21 
For comparison, in Fei. 76 Ge the ordered moments lie in 
the a&-plane, and an average moment of 1.56±0.2/is/Fe 
was reported based on neutron diffraction. 27 There have 
been contradictory reports regarding a variation of the 
moments between the two Fe sites in Fe 2 - 2 ;Ge. 27-29 The 
more recent Mossbauer results have suggested that a 
larger moment resides on the Fe(2) site, and perhaps 
similar experiments on Fe 3 GeTe 2 may provide additional 
insight into the roles of the different Fe environments. 
The difference in the easy axis between Fe 3 GeTe 2 and 
Fei. 67 Ge is most likely caused by the increased chemi¬ 
cal anisotropy associated with the inclusion of the Te-Te 
layer into the FesGeTe 2 structure. 

We confirmed the ferromagnetic ordering temperature 
and orientation of the moments by tracking the 010 and 
002 Bragg peaks, as shown in Fig. 6. The power law fit 
between 80 K and 180 K in Fig. 6 yields T c = 148(3) K, 
which is consistent with the bulk magnetization measure¬ 
ments. The intensity of the 010 peak increases when 
the moments lie perpendicular to the 010 scattering vec¬ 
tor. Therefore, the increase in intensity below ss 150 K 
shown in Fig. 6(a) demonstrates that the moments do 


not lie along the 6-axis, and when combined with the 
temperature-independent behavior of 002 intensity we 
verify that the moments lie along c. This shows that 
there is not any significant spin canting or reorientation 
as temperature decreases, which is also confirmed with 
the powder neutron diffraction measurements. In addi¬ 
tion, we did not observe any significant change in the 
nuclear structure across the magnetic transition. 
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FIG. 7. (color online). Neutron powder diffraction at (a) 
200 K and (b) 60 K with refinements included and samples 
labeled by the refined compositions. The reflections from top 
to bottom correspond to the structural, aluminum can and 
magnetic reflections. The most intense magnetic signal was 
observed in the region shown in figures (c) and (d). Note the 
much increased magnetic intensity in the Fe 2 . 9 GeTe 2 sample 
compared to the more Fe-deficient Fe 2 . 7 iGeTe 2 - 


During the single crystal neutron diffraction measure¬ 
ments, we scanned along H and L and found no addi¬ 
tional non-integer (H,K,L) intensity through the mag¬ 
netic transition of 150 K. Additionally, in the neutron 
powder diffraction data shown in Fig. 7, there is no evi¬ 
dence of scattering at non-integer (HKL) positions in go¬ 
ing from above the Curie temperature Tc to below Tc- 
This strongly suggests that there is no deviation from the 
observed ferromagnetic structure for both Fe-rich and Fe- 
deficient samples. We did notice some diffuse scattering 
beneath the Bragg peaks in our single crystal neutron 
diffraction, and this scattering was persistent at temper¬ 
atures well above the magnetic ordering. As such, this is 
likely related to crystalline defects and requires further 
investigation (note that stacking faults were not observed 
in our STEM). 



















TABLE III. Summary of refinements of powder neutron diffraction data. Refined compositions are the result of refining the 
occupation of the Fe(2) position. Refined moments are aligned along the c-axis. 


Nominal Composition 

T 

(K) 

Refined Fe 

a 

(A) 

c 

(A) 

moment Fe(l) moment Fe(2) Rp 
(Ms) (/^b) 

Fe 3 GeTe 2 

250 K 

2.904(8) 

4.01749(8) 16.33990(8) 

- 

- 

3.02 

FesGeTe2 

60 K 

2.904(8) 

4.00938(3) 

16.2850(2) 

2.18(10) 

1.54(10) 

3.41 

Fe2.75GeTe2 

2001< 

2.71(3) 

3.95001(7) 

16.4019(4) 

- 

- 

3.42 

Fe2.75GeTe2 

60 K 

2.71(3) 

3.93628(4) 

16.3535(3) 

1.4(1) 

1.4(1) 

5.25 


Neutron powder diffraction was performed on poly- 
crystalline samples of nominal compositions Fe 3 GeTe 2 
and Fe 2 . 75 GeTe 2 (see Fig. 7 where refined compositions 
are utilized). In both cases, no deviation from the 
magnetic structure obtained from single crystal neu¬ 
tron diffraction was observed, indicting similar ferromag¬ 
netism regardless of Fe concentration. The refinement 
results are summarized in Table III. The Fe-rich sample, 
with a refined composition of Fe 2 .goGeTe 2 , was found to 
have a stronger magnetic contribution to the Bragg peaks 
and the refinement yielded a larger ordered moment on 
Fe(l) than on Fe(2) (the raw data for this sample are 
shown in the Supplemental Materials). Fixing the mo¬ 
ments to be the same on the Fe(l) and Fe(2) sites did 
not provide suitable fits to the data for this sample. This 
result is consistent with a recent publication containing 
neutron powder diffraction on a sample of Fe 2 . 9 oGeTe 2 - 21 
The behavior is different in the Fe-deficient sample, which 
has a refined composition of Fe 2 . 7 iGeo. 95 Te 2 . In this sam¬ 
ple, the magnetic contribution to the diffraction peaks is 
smaller and appears similar for both Fe(l) and Fe(2). 
The refined moment is slightly larger than that obtained 
on the single crystal (similar composition), but the trend 
for similar moments on the Fe(l) and Fe(2) sites is con¬ 
sistent. This suggests that vacancies on the Fe(2) site 
suppress the magnetism on both Fe sites. This could be 
related to the changes in structural parameters, though 
disorder and/or dilution effects may also be dominant. 


C. Magnetization 

Magnetization measurements were performed on the 
polycrystalline materials to determine their Curie tem¬ 
peratures and correlate this with structure and compo¬ 
sition. Results of the temperature-dependent magneti¬ 
zation M measurements are shown in Fig. 8(a), where 
refined Fe contents are used in the legend. We have used 
Gaussian-CGS units for the magnetization results. 

The data in Fig. 8(a) demonstrate that the Curie tem¬ 
perature decreases with decreasing Fe content. The re¬ 
duced Xc with increasing vacancies on the Fe(2) site may 
be caused by a disruption of the magnetic exchange with 
increasing disorder and magnetic dilution via vacancies. 
There may also be a structural component, as we observe 
that the reduced Tq correlates with the expansion of c, 
an increase in Fe(l)-Fe(l) bond distance, and decrease in 
Fe(l)-Fe(2) bond distances. A similar reduction in Tq is 



T( K) 


FIG. 8. (color online) (a) Temperature dependence of the 
magnetization M divided by applied field H for polycrys¬ 
talline Fe 3 - a ,GeTe 2 materials with various Fe concentrations 
as indicated by the refined Fe content in the legend, and 
(b) anisotropic magnetization of flux-grown Fe2.76Geo.g4Te2. 
All data were collected upon cooling in an applied field of 
77=100 Oe. 


observed for Fe 2 - K Ge materials as x increases, though, 
where a more typical decrease in the lattice parameters is 
observed with increasing x. 19 Investigating the pressure- 
dependence of Tc or the anisotropy of the magnetic ex¬ 
citation spectra may provide further insight, as would 
theoretical calculations into the dependence of T c on c. 

Anisotropic magnetization data for single crystalline 
Fe 3 _a,GeTe 2 are shown in Fig. 8(b). These data demon¬ 
strate that the easy axis for magnetization is along the 
crystallographic c-axis, which is consistent with our neu¬ 
tron diffraction and Ref. 1. The shape of M(T) evolves 
with decreasing Fe content. In the Fe-rich samples, M{T) 
increases smoothly while cooling. In contrast, for the 
most Fe-deficient sample, M(T) has a kink-containing 
shape near Tq and essentially plateaus slightly below Tc- 
The M(T ) data reported on the vapor transport crystals, 
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with Tc ~ 220 K, has a temperature dependence similar 
to that shown in Fig. 8(b). Therefore, the non-power law 
behavior observed near Tq in the single crystals may be 
linked to domain wall formation and movement. It is cer¬ 
tainly possible that different growth conditions produce 
different Fe(2)/Fe(3) contents, which could also modify 
the properties independently. We also confirmed that 
Ni 3 _ a; GeTe 2 appears to be a Pauli paramagnet. 



H (kOe) 

FIG. 9. (color online), (a) Magnetization loops at 2K for 
three compositions of polycrystalline Fe 3 _ x GeTe 2 (refined 
values provided), and the inset shows the magnetization at 
2 K and 60 kOe for polycrystalline samples, (b) Magnetiza¬ 
tion versus applied field for single crystalline Fe 2 . 76 Geo. 94 Te 2 
at various temperatures, with the anisotropy demonstrated at 
2 K. 

Isothermal magnetization data are shown in Fig. 9. 
In panel (a), data for three polycrystalline samples are 
shown while panel (b) contains data for an oriented sin¬ 
gle crystal. The isothermal magnetization measurements 
for polycrystalline samples reveal some interesting trends 
with composition. It is clear that the Fe-deficient samples 
have a much lower remanent magnetization and very lit¬ 
tle coercivity. Also, the moment essentially saturates at 
high fields for the Fe-deficient samples while a linear rise 
with field is observed for the Fe-rich samples. The linear 
increase at higher fields for Fe 2 . 97 GeTe 2 may be related 
to an increased anisotropy field at higher Fe content. Our 


a (A) 

3.94 3.96 3.98 4.00 4.02 4.04 



c(A) 


FIG. 10. (color online). Curie temperature as a function of 
lattice parameters for Fe 3 _:,;GeTe 2 . 


low-Fe content crystals have an estimated anisotropy field 
of 15kOe (or 1.5 T), as shown in Fig. 9(b). An anisotropy 
field of «5 T was observed at 10 K for crystals grown by 
vapor transport, 1 and the reported Tq ~ 220 K suggests 
a large Fe content in those crystals. The linear rise at 
large fields may be also be from some paramagnetic ions, 
such as an impurity or interlayer Fe. We note that our 
Fe 2 . 76 Geo. 94 Te 2 crystals did not reveal any unexpected 
behavior when measured to 120 kOe. 

The ‘saturation magnetization’ of all polycrystalline 
Fe 3 _a,GeTe 2 samples is shown as an inset in Fig. 9(a) 
and values are listed in the summary of samples provided 
as Table IV. The saturation magnetization is taken as 
the value of the magnetization (in fig/ Fe) obtained at 
2 K and 60 kOe (refined compositions are used for unit 
conversion). This does not represent a true saturation 
magnetization for all samples due to the linear increase 
in M at large H. We clearly see that higher Tq (higher 
Fe content) correlates with larger induced moments for 
a given T at large applied fields, consistent with our 
powder neutron diffraction results. T c was defined us¬ 
ing the intercept of the steepest tangent. The effective 
moments calculated from the susceptibility (x) between 
250 and 360 K were found to vary between 3.9(2) and 
4.9(l)p,s/Fe, with higher T c generally corresponding to 
larger effective moments. We used a standard Curie- 
Weiss law (x = C/(T — Tew)) and fit data collected 
on cooling in an applied field of 1 kOe. Curie-Weiss tem¬ 
peratures Tew obtained from these fits agreed well with 
the Curie temperatures obtained from measurements at 
lower fields. 

We calculated the Rhodes-Wohlfarth ratio (RWR) 
for our polycrystalline samples. The RWR provides a 
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TABLE IV. Summary of samples and magnetic properties for polycrystalline Fe 3 _ a; GeTe 2 . Rietveld compositions originate in 
refinement of powder diffraction data while EDS values are for large area EDS analysis that may include impurities. 

Fe composition a c T c p, sa t Me// Tew 

nominal Rietveld EDS (A) (A) (K) (mb/ Fe) (mb/ Fe) (K) 


3.10 

2.97(2) 

3.07(2) 4.0244(1) 16.3293(5) 232 

1.31 

4.9(1) 

225.6(2) 

3.00 

2.96(1) 

3.06(3) 4.0184(1) 16.3344(4) 226 

1.32 

4.8(1) 

221.7(1) 

2.90 

2.92(1) 

2.92(3) 4.0013(1) 16.3494(4) 220 

1.28 

4.4(1) 

218.9(1) 

2.85 

2.87(1) 

2.88(2) 3.9840(1) 16.3676(3) 203 

1.27 

4.4(1) 

203.2(1) 

2.80 

2.78(1) 

2.81(2) 3.9564(1) 16.3947(4) 154 

1.08 

4.5(1) 

143.5(3) 

2.75 

2.69(4) 

2.79(3) 3.9475(4) 16.405(1) 140 

1.04 

3.9(2) 

142.6(2) 


quick means to characterize the degree to which a mag¬ 
netic moment is localized. RWR is defined as RWR 
= Pc/Ps , with p c obtained from the effective moment 
Pc(p c + 2) = p^. 30,31 Physically, p c is the saturation 
moment expected from the effective moment calculated 
from the susceptibility in the paramagnetic phase (as¬ 
suming Curie-Weiss behavior) and p s is the saturation 
moment obtained in the ordered state. RWR = 1 for 
localized systems and is larger in an itinerant system, 
with the ratio increasing as Tq decreases for itinerant 
systems. Here, we take p s as the magnetization obtained 
at 2 K and 60 kOe and calculate RWR values between 2.7 
and 3.4 for our polycrystalline samples. These values are 
fairly similar to the RWR = 3.8 reported in Ref. 1. 

The observation of RWR > 1 in compounds with a low 
Curie temperature (Tc < 500 K) suggests itinerant fer¬ 
romagnetism is likely present. In comparison to the val¬ 
ues tabulated by Wohlfarth in 1978 and Moriya in 1979, 
these Fe 3 _a;GeTe 2 compounds lie in the region between 
localized and itinerant ferromagnetism. 30,31 We do not 
observe a strong magnetoelastic effect at the transition, 
which is expected for large RWR (itinerant) systems. 30 
While we do not observe a strong increase in the RWR 
values as T c decreases, this can likely be attributed to 
the influence of vacancies on the structure and mag¬ 
netism (multiple effects influencing both Tq and RWR). 
Future measurements of the spin-waves via inelastic neu¬ 
tron scattering will aid in addressing the itinerant nature 
of this system, as would theoretical or experimental stud¬ 
ies into the influence of pressure on Tq. 


(a) 
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zL 



V 0 H(T) 



FIG. 11. (color online) (a) Hall resistance as a function of ap¬ 
plied field suggests p-type conduction in Fe2.76Geo.94Te2 and 
Ni 3 _ x GeTe 2 crystals, while in (b) the negative Seebeck coeffi¬ 
cients indicates n-type conduction. The dashed line in panel 
(b) indicates the Curie temperature of the Fe2.76Geo.94Te2 
crystal. 


We have summarized the magnetization data as a plot 
correlating the lattice parameters with the Curie tem¬ 
perature, shown in Fig. 10. In addition, Table IV pro¬ 
vides a summary of sample compositions determined us¬ 
ing various methods as well as the magnetic properties. 
Tc clearly decreases as Fe vacancies are introduced and 
the lattice responds with a decrease in the in-plane lat¬ 
tice parameter and a slight expansion along c. These 
results can be used as a guide to predict the composition 
needed to obtain a particular Tq or as a means to expedite 
characterization, and must be considered when perform¬ 
ing measurements on ultra-thin samples. Figure 10 shows 
that the flux-grown crystals do not behave unexpectedly 
based on the behavior of polycrystalline materials. 


D. Hall Effect and Seebeck Coefficients 

We have utilized Hall effect and thermoelectric mea¬ 
surements to characterize the in-plane electrical trans¬ 
port in Fe2.76Geo.94Te2 and Ni 3 _ x GeTe 2 crystals. These 
results, summarized in Fig. 11, reveal that both systems 
likely have multiple carrier types contributing to conduc¬ 
tion. We observe positive Hall coefficients and negative 
Seebeck coefficients, which would independently suggest 
p-type and n-type conduction, respectively. The linear 
dependence of the Hall resistance on applied field pro¬ 
hibits a detailed analysis aimed at determining the con- 
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tribution of each band/carrier-type. Theoretical calcu¬ 
lations would provide additional insight into the origin 
of this apparent multi-carrier transport. We note that 
our flux-grown crystals of Ni 3 _ a; GeTe 2 likely have a^O.O 
based on EDS measurements. 

At room temperature, the Hall coefficient f?H of 
Ni 3 _ a ,GeTe 2 is about an order of magnitude smaller than 
that of Fe 2 . 76 Geo. 94 Te 2 , which translates to a larger hole 
concentration in Ni 3 _a,GeTe 2 if a single-carrier model is 
used. Specifically, at 300 K the Hall carrier density ?4 h = 
l/i?ne is approximately 1.8xl0 22 cm -3 for Ni 3 _ a; GeTe 2 
and approximately 1.9xl0 21 cm -3 for Fe 2 . 76 Geo. 94 Te 2 - If 
both holes and electrons are present, as suggested by 
these results, the Hall coefficients could be artificially 
reduced and the carrier concentrations reported would 
be upper-limits to the actual number of holes in the sys¬ 
tem. A more complete compensation of charge carri¬ 
ers may be responsible for the smaller Hall coefficient 
of Ni 3 _ a; GeTe 2 , though we found Ni 3 _ ;c GeTe 2 to have 
about an order of magnitude lower electrical resistivity 
(Supplemental Materials). As discussed in the Supple¬ 
mental Materials, the Hall coefficient of Fe 2 . 76 Geo. 94 Te 2 
is strongly influenced by an anomalous Hall contribution 
below «200K. Thermal conductivity and specific heat 
data for Fe2.76Geo.94Te2 and Ni 3 _ a ,GeTe 2 are also pre¬ 
sented in the Supplemental Materials. 32 

IV. SUMMARY 

The availability of Fe 3 _a,GeTe 2 and other recently- 
developed, cleavable ferromagnets provides a starting 
point for the development of magnetically-active van 
der Waals heterostructures. While such architectures 


will likely be designed to investigate specific physics or 
functionality, fundamental investigations will almost cer¬ 
tainly yield unique spin structures or magnetotransport. 
This work has demonstrated that the itinerant ferromag¬ 
netism in Fe 3 _ a; GeTe 2 can be tuned by controlling the 
Fe content. All manifestations of the magnetic interac¬ 
tions, from the Curie temperature to the local moment, 
are reduced as Fe vacancies are created and the in-plane 
lattice parameter decreases. By mapping the magnetic 
phase diagram of Fe 3 _ a; GeTe 2 , this work has provided a 
foundation for future studies examining the influence of 
dimensionality on the magnetism in these van der Waals 
bonded materials. While we have shown that the mag¬ 
netic behavior can be controlled through total Fe content, 
chemical substitutions or intercalation may provide ad¬ 
ditional control over the magnetism and physical proper¬ 
ties of this layered material. Future work in this area will 
need to examine the stability of ultra-thin Fe 3 _a,GeTe 2 as 
a function of Fe content. In addition, experiments under 
pressure and theoretical calculations will likely provide 
valuable information regarding the interactions between 
the electronic and magnetic structures of these multi¬ 
carrier, itinerant ferromagnets. 
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A. Supplemental Materials 

Neutron powder diffraction data collected on a sam¬ 
ple of nominal composition Fe3GeTe2 are shown in 
Fig. SI. For this sample, the refined composition was 
Fe 2 .904(g)GeTe2. Refinement results are summarized in 
the main text. 

With the production of large single crystals, we were 
able to perform in-plane thermal and thermoelectric 
transport measurements. Figure S 2 presents the elec¬ 
trical resistivity, Seebeck coefficient, and thermal con¬ 
ductivity of our flux-grown Fe 2 .76(4) Ge 0 .94(4)Te 2 crystals 
(in-plane transport). This composition is obtained from 
refinement of single crystal neutron diffraction data. To 
facilitate a comparison, we have used the composition 
Ni2.4o(4)Ge! 0 i(3)Te 2 .oo(3) obtained from EDS. The elec¬ 
trical resistivity is about an order of magnitude lower 
for Ni 2 .4oGeTe2 than for Fe 2 . 76 Geo. 94 Te 2 (see caption of 
Fig. S 2 ). The Ni 2 .4oGeTe 2 sample also has a slightly 
lower residual resistivity ratio (RRR), as demonstrated 
by the plot of p(T)/p( 300 K) in Fig. S 2 (a). Transition 
metal vacancies and/or associated displacements of Ge 
likely provide significant charge carrier scattering that 
leads to small RRR in both of these systems. 
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FIG. SI. (color online) Neutron powder diffraction at (a) 
250 K and (b) 60 K with refinements included and sample la¬ 
beled by the refined composition. The reflections from top 
to bottom correspond to the structural, Aluminum can and 
magnetic reflections. 


The Seebeck coefficient (a) is small and negative for 
both samples. The negative value implies electrons dom¬ 
inate conduction, and the smaller value for Ni 2 .4oGeTe2 
would imply a higher concentration of charge carriers 
(consistent with lower p) or a more complete compensa¬ 
tion of electrons/holes. Based on the Hall data discussed 
below, these appear to be multi-carrier metals. 

The estimated lattice thermal conductivity Ki a t is sim¬ 
ilar for Fe 2 . 76 Geo. 94 Te 2 and Ni 2 . 4 oGeTe 2 , as shown in 
Fig. S 2 (c). These Ki at values were obtained using the 
Wiedemann-Franz law to estimate an electronic contri¬ 
bution K e to the total thermal conductivity k\ the degen¬ 
erate limit of the Lorenz number was assumed. The val¬ 
ues of Ki at are similar across the entire temperature range 
investigated and a low T maximum is not observed. The 
temperature dependence suggests vacancies likely domi¬ 
nate phonon scattering rates, and despite apparently dif¬ 
ferent transition metal contents (based on our EDS re¬ 
sults) the net result is a similar ni at . It is also possible 
that phonons are scattered by charge carriers at low T. 

The electrical properties respond to the ferromagnetic 
ordering in Fe2.76Geo.94Te2, as observed in Fig. S 2 (a,b) 
where the Curie temperature is indicated by the dashed 
line. Below Tq, the electrical resistivity and Seebeck co¬ 
efficient also begins to decrease more rapidly. A decrease 
in p below Tc is commonly understood as a reduction in 
spin disorder scattering when the moments order. The 
source for a decrease in the Seebeck coefficient is less 
clear. The Seebeck coefficient is influenced by the scatter¬ 
ing mechanisms as well as the carrier concentration and 
shape of the Fermi surface, the latter of which may be 
influenced by the magnetic ordering. The magnetic scat¬ 
tering may asymmetrically influence the contributions of 
the various charge carriers, thereby influencing the See¬ 
beck coefficient. We have not noticed a strong response 
of the lattice to the magnetic ordering in Fe 2 .76Geo.94Te 2 . 
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FIG. S2. (color online) In-plane electrical and thermal trans¬ 
port for Fe2.76Geo.94Te2 and Ni2.4oGeTe2 crystals, (a) Nor¬ 
malized electrical resistivity, (b) Seebeck coefficient and (c) 
thermal conductivity with estimates for the lattice and elec¬ 
tronic contributions shown. The closed black symbols are for 
Fe 2 . 76 Geo. 94 Te 2 (p(300K)=3.0xl0 _4 f2-cm) and the open blue 
symbols represent Ni 2 .4oGeTe 2 (p(300K)=9.50xl0“ s n-cm). 


We performed Hall effect measurements to further 
characterize the electrical behavior of these systems. The 
Hall coefficient is positive for both Fe 2 . 76 Geo. 94 Te 2 and 
Ni 2 . 4 oGeTe 2 , and the Hall resistance pn is linear with 
magnetic fields up to at least 8 T at 300 K (Fig. S3(a)). At 
room temperature, the Hall coefficient Ru of Ni 2 . 4 oGeTe 2 
is about an order of magnitude smaller than that of 
Fe 2 . 76 Geo. 94 Te 2 , which translates to a larger carrier con¬ 
centration in Ni 2 . 4 oGeTe 2 if a single-carrier model is used. 
Specifically, at 300 K the Hall carrier density «h = 1/Rtie 
is ssl.8x 10 22 cm -3 for Ni 2 . 4 oGeTe 2 and «1.4x 10 22 cm -3 
for Fe 2 . 76 Geo. 94 Te 2 using the linear fits shown in Fig¬ 
ure S3(a). 

The positive Hall coefficient suggests the dominant car¬ 
riers are holes while the negative sign of the Seebeck co¬ 
efficient suggests the dominant charge carriers are elec¬ 
trons. A detailed analysis of the contributions of each 
band is prohibited, however, due to the linearity of the 


FIG. S3. Hall effect data at (a) 300 K for Fe2.76Geo.94Te2 
and Ni 2 . 4 oGeTe 2 and in (b,c) the influence of an anomalous 
component brought about by the ferromagnetic ordering in 
Fe2.76Geo.94Te2 is demonstrated. 


Hall resistance with magnetic field. If both holes and 
electrons are present, as suggested by these results, the 
Hall coefficients could be artificially reduced and the car¬ 
rier concentrations reported would be upper-limits to the 
actual number of holes in the system. Ni 2 . 4 oGeTe 2 is 
more conductive, despite apparently having more defects, 
and thus an increase in the absolute number of charge 
carriers relative to Fe 2 . 76 Geo. 94 Te 2 is likely. 

The Hall effect of Fe 2 . 76 Geo. 94 Te 2 is strongly influ¬ 
enced by the anomalous Hall contribution, as shown in 
Fig. S3(b,c). For our crystals, with a Curie temperature 
of ss 150 K, an influence of the anomalous Hall effect is 
observed below approximately 200 K. This is due to the 
strong polarization of the Fe moments with increasing 
field and decreasing temperature. The Hall data follow 
the field dependence of the magnetization, which demon¬ 
strates that the non-linearity of pn is not due to multi¬ 
band effects. The regular and anomalous Hall coefficients 
have the same sign, and the current data are insufficient 
to analyze in detail due to the small contribution from 
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FIG. S4. (color online). Specific heat capacity of 
Fe2.76Geo.94Te2 and Ni2.4oGeTe2, with inset showing low T 
behavior and fits (thin solid lines). A small anomaly is ob¬ 
served in Fe2.76Geo.94Te2 near Tc«150K. 


the regular Hall coefficient as well as the non-saturating 
magnetization at high fields. Qualitatively different re¬ 
sults were obtained for the regular Hall coefficient when 
a detailed analysis was performed on data collected for 
different crystals, which had similar room temperature 
Hall coefficients (likely due to minor variations in mag¬ 
netization between the crystals). 


Specific heat measurements were performed on sin¬ 
gle crystals of Fe2.76Geo.94Te2 and Ni 2 .4oGeTe2, and 
the results are shown in Fig. S4. A small anomaly is 
present in the region of the ferromagnetic transition of 
Fe2.76Geo.94Te2 while data for Ni 2 .4oGeTe2 are smooth 
across the entire temperature range. The measured 
values are approaching the high temperature limit of 
3fcs/atom at 220 K for both materials. This is consistent 
with the Debye temperatures obtained from the low T 
data, which were 0u=224K and 234 K for Ni 2 . 4 oGeTe 2 
and Fe 2 . 76 Ge 0 . 94 Te 2 , respectively. The slightly smaller 
0 d for Ni 2 . 4 oGeTe 2 could be due to increased softening 
associated with a higher carrier density or a higher va¬ 
cancy concentration, though potential error associated 
with the sample compositions precludes such a conclu¬ 
sion. At low temperatures, the Debye temperature Op 
is obtained from a plot of Cp/T versus T 2 plot where 
the slope is 12 ^ e y i JV °* and N a t is the number of atoms 
per formula unit. 

The electronic coefficient to the specific heat 
is significantly larger for Fe 2 . 76 Geo. 94 Te 2 than for 
Ni 2 . 4 oGeTe 2 . We obtain a Sommerfeld coeffi¬ 
cient 7 = 132.8 mJ/mol/K 2 for Fe 2 . 76 Ge 0 . 94 Te 2 and 
10.2 mJ/mol/K 2 for Ni 2 . , 4 oGeTe 2 . The large 7 for 
Fe 2 . 76 Geo. 94 Te 2 is either due to a mass enhancement from 
correlations or due to a contribution from spin fluctua¬ 
tions. We note that our 7 is similar to that reported for 
vapor transport grown Fe 3 GeTe 2 crystals . 81 
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